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Abstract 
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I. INTRODUCTION 


Interest in spin physics has grown in the last few years thanks to the various polarized DIS 
experiments and to existing proposals for future polarized beam experiments. One proposal 
includes the construction of a polarized ep collider at DESY (HERA-Y) running at an 
energy of ~ 40 GeV, which will complement and extend the results to be obtained at the 
Relativistic Heavy Ion Collider (RHIC) at Brookhaven. Polarized deep-inelastic scattering 
experiments have shed light on the relative contributions of the nucleon constituents to their 
overall spin. However, the largest theoretical uncertainties are from: (1) the polarized gluon 
size and shape as a function of x being unknown and (2) the relative size of the polarized 
sea distributions being uncertain. |]l[ 

In a recent paper, we proposed three new parametrizations for the spin dependent 
parton distributions for the proton. Our distributions, which are available in both leading 
order (LO) and next-to-leading order (NLO), were evolved directly in x-space. In ref. 
we presented details of the models used to obtain the input distributions and compared 
the evolved structure functions with the available data. All currently available data on the 
nucleon spin structure come from deep inelastic scattering experiments, and therefore do not 
contain direct information on individual parton distributions such as the gluon, AG{x, 
or strange sea, As{x,Q‘^), for example. It is expected that such detailed information will 
come from other experiments such as those proposed at RHIC and HERA-Y. 

As a pp collider, RHIC will be endowed with a very high integrated luminosity (up to 800 
pb“^ at 500 GeV) and will be spanning a center of mass energy range between 50 to 500 GeV. 
One of the main programs at RHIC will be to determine the size and shape of the polarized 
parton distributions which, at the moment, suffer from signihcant model dependence. There 
are large uncertainties in both the gluon contribution, AG, and the sea quark distributions. 

In this paper we continue the work of ref. |0] by providing details of the x-space evolu¬ 
tion of the parton distributions. We also compare our evolved distributions with another 
parametrization in order to highlight the similarities and differences between them. Finally, 
we use our evolved parton distributions to predict cross sections and asymmetries for direct 
photon and single jet production at RHIC and HERA-Y energies. 


II. INPUT DISTRIBUTIONS 

In reference p|, we constructed three sets of polarized parton distributions for the va¬ 
lence and all light sea flavors of quarks. Each set corresponds to a particular model of 
the polarized gluons, motivated by different physical assumptions. All of the x-dependent 
polarized quark distributions are generated from the unpolarized distributions and some 
basic theoretical assumptions. The valence up and down distributions are generated from 
a broken SU(6) model of the polarized valence quarks. For the sea distributions, we built 
in positivity constraints and incorporated the integrated and x-dependent polarized deep- 
inelastic scattering data. The three models of the polarized gluons used to generate the 
distributions are assumed small at the relatively low values of the data. One model 
gives a positively polarized gluon, generated from the unpolarized glue: (AG = xG), one 
is an unpolarized gluon (AG = 0), and the third is based upon an instanton model, which 
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integrates to a slightly negative value. All three sets of distributions are in good agreement 
with the data, an indication that the polarized distributions are still not well determined. 

All of our distributions are generated at = 1.0 GeV^ and are evolved up to the 
values necessary for predicting the appropriate spin observables. Evolution for the singlet 
and non-singlet functions is done in x-space to the NLO level. We have assumed three 
quark flavors up to the appropriate value for charm production, whereby we smoothly 
introduce the fourth flavor to the evolution equations. Details will be discussed briefly in 
the next section. 

The distributions we use here are strongly embedded in sound theoretical assumptions, 
while giving good agreement with all of the latest polarized DIS data. In order to determine 
the nature of the polarized distributions to a hner degree of detail, we are using them to 
generate predictions of the spin asymmetries for direct photon and jet production. These 
processes are feasible for the RHIC and HERA experimental groups and there is general 
agreement that they are pivotal in helping to determine the nature of the polarized gluon 
distribution. 


III. A-SPACE EVOLUTION IN LEADING AND NEXT-TO-LEADING ORDER 


Given the polarized parton distributions at some reference scale Ql, they can be predicted 
at a higher scale by solving the polarized DGLAP evolution equations. For the non-singlet 
sector the evolution equation is given by 


dAqNs{x,Q^) as(Q‘^ 


d\nQ^ 
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-AP^s{x) * AqNs{x, Q^), 


(3.1) 


where AqNsix,Q‘^) is the polarized non-singlet parton distribution function, AP^g{x) are 
polarized non-singlet splitting functions, asiQ"^) is the strong coupling constant and * rep¬ 
resents the convolution 


fix) * gix) = j —/ ("-"j giy). 
Jx V \v 


y \y, 

For the polarized singlet sector, with the quark singlet dehned by 

AS = 


the evolution equations are 
dAS(a;,g2) 


d In 271 

dAGix,Q^) asiQ^) 


APg,{x) * AS(a;, Q^) + APqg{x) * AG(a;, Q^) 
APggix) * AEiX, g2) + APggix) * AG (X, Q^] 


dlng^ 271 

The splitting functions have the perturbative expansion 


AP,gix) = AP}^\x) + 


CXs i^Q 


Ai’«(x) + ... 




(3.2) 


(3.3) 


(3.4) 


(3.5) 
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The first term is the leading order (LO) contribution and the second term the next-to- 
leading order (NLO) contribution to the splitting kernels, which were calculated in ref. [Q. 
In leading order we use the expression 




dvr 

/5oln(QVA2)’ 


(3.6) 


and in NLO 




dvr 

/?oln(QVA2) 


A lnln(gVA2) 


(3.7) 


for the strong coupling constant. The constants /do and Pi are dehned by 

^-N,Nf - 2CFNf, 


Po - yA^c - 2 A^/, 




(3.8) 


with the color constants N^. = 3 and Cp = d/3, and Nf is the number of active quark flavors. 

In order to solve the DGLAP evolution equations directly in a:-space we use the ansatz 
proposed in ref. 0]. We assume a solution of the form 


Aq{x, Q" 


^ ^ Anjx) 


n=0 


n\ 


'apQ^r 

MQl). 


+ as(g^) X! 

n=0 






In" 


n\ 


aJ,Ql 


(3.9) 


where Aq = {Aq ns, AH, AG). Since the sum starts at n = 0, boundary conditions, 
q{x,Ql) = qo{x), can be incorporated into the equations. When Eq. (|3.9|) is inserted 
into the DGLAP equations and all third-order terms are dropped, the relations 


Ao(x) = qQ{x) 

Bo{x) = 0 

An+i{x) = * An{x) 

Po 

Bn+i{x) = (-5(1 -x) - ^P^°\x)) * Bn{x) - ^P^^\x) * An{x) - -J^An+px) (3.10) 

Po 'xPo dTrPo 

are obtained for the An{x) and Bn{x). The functions P^^’^px) are the LO and NLO parts 
of the splitting functions, while the terms An(x) and Bn(x) represent the LO and NLO 
coefficients respectively. If a LO evolution is desired then only the An(x) are calculated and 
the LO expression for is used. These equations can be solved numerically for any value 
of X, and the coefficients substituted into Eq. ( |3.9|) to obtain the distributions Aq{x, Qp at 
any value of Q^. Details of the input distributions and expressions for the coefficients are 
given in the Appendix. 


IV. LARGE Pt processes AT RHIC AND HERA-N 
A. Experimental Parameters 

The polarization experiments planned for RHIG show great promise for the high energy 
community. With polarized proton beams of 70% polarization and a planned luminosity of 
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C = 2 X 10^^/cm“^/sec“^, both direct -7 production and jet production can be done in a 
kinematic region where we can gain tremendous insight into the polarized gluon distribution. 
These are part of the planned set of experiments for both the STAR and PHENIX detectors 
at RHIC. 0 If the experiments are able to obtain an integrated luminosity, presently as¬ 
sumed to be 320 pb~^ at y/s = 200 GeV and 800 pb~^ at ^/s = 500 GeV, the resulting data 
should be enough to distinguish between various polarized gluon models which have been 
proposed, 

The STAR detector will have a wide angular coverage to cover a large kinematic region 
for Jet production. Their range of x will extend from about 0.02 < x < 0.30 (±0.02) with 
an approximate uncertainty of 6{AG) ~ ±0.05. 

The PHENIX detector will measure high pt prompt 7 ’s in the kinematic range of 10 < 
Pt < 30 GeV at y/s = 200 GeV {x < 0.25) and 10 < pr < 40 GeV at y/s = 500 GeV 
{x < 0.15). The angular coverage is not as large as STAR but it has hner granularity. 
The estimated uncertainty is 5{^} ~ ±0.05 —> 0.30, increasing linearly with px- The 
combination of the two detectors will be valuable in covering a wide kinematic range with 
reasonable error bars for cross sections and asymmetries. 

Tests are presently being performed to polarize the proton beam at HERA. The success 
of these tests will provide a unique kinematic region for measuring asymmetries in direct 
7 and Jet production. [§] These experiments will be performed at y/s = 40 GeV, giving a 
good sensitivity to distinguishing between gluon model predictions for the corresponding 
asymmetries. This will be discussed later. 

All of the above mentioned experimental groups plan to measure the observables in both 
direct -7 and Jet production. Present plans call for the RHIG experiments to be performed in 
the next hve years, while the HERA experiments would hopefully follow shortly thereafter. 
Both sets of experiments are important, as the lower center-of-mass energies appear to be 
somewhat better for distinguishing the size of the polarized gluons. Thus, the combination 
of a wide kinematic range and reasonable uncertainties should give a good indication of the 
size and shape of the polarized gluon distribution. 


B. Theoretical Background 


1. Prompt Photon Production 


The prompt photon cross section promises to be one of the most useful for measuring 
AG at RHIG and HERA-V since it is dominated by initial qg scattering. Gontributions 
to the prompt photon cross section are usually separated into two classes in both LO 0 
and NLO. |T0HT2| There are the so-called direct processes, ab —>■ 7 c in LO and ab —>■ 'ycd, 
in NLO, a, b, c and d referring to partons, where the photon is produced directly in the 
hard scattering. In addition there are the fragmentation contributions where the photon 
is produced via bremsstrahlung off a hnal state quark or gluon, ab —>• cd{e) followed by 
c —>■ 7 ± A for instance. 

In LO, 0{aas), the direct subprocesses contributing to the cross section are 


Q9 

qq 


7Q 

19- 


(4.1) 
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In addition there are the fragmentation processes 


qg qg 

qq qq 
qq' qq' 
qq^ qq 

qg qg 
qq gg 
gg gg 

gg qq (4.2) 


where one of the hnal state partons fragments to produce the photon, he., q{g) 
In the direct processes in LO, the differential cross section is given by 

„ dAai° 1 _ r'' dv 


'1 + X. 






vw 1 


Ivw/v V dv 


5(1 — w) (4.3) 


where S = {Pi + P 2 )^ V" = 1 + T/S, W = -U/{T + S),v = l + i/s, w = -u/{i + s), 
s = X 1 X 2 S, and T = (Pi — P^)^ and U = (P 2 — Py)^. As usual the Mandelstam variables 
are dehned in the upper case for the hadron-hadron system and in lower case in the parton- 
parton system. Pi and P 2 are the momenta of the incoming hadrons and Afl{xi, M^) and 
Af 2 {x 2 , M^) represent the respective probabilities of hnding parton i and j in hadrons 1 
and 2 with momentum fractions Xi and X 2 at scale M^. The quantity dA&ij^^/dv is the 
hard subprocess cross section which is calculable in perturbative QCD. 

For the fragmentation processes, the differential cross section is given by 




d^P-y 




1 dAa, 

V dv 


dz dv 

VWjz 1 — V 


ttS f^^Ji-v+vw 

^^^'6{l-w)Dy/i{z,Mf), 


—Afi{xi,M^)Afi{x2,M^) 

Wwjvz W 


(4.4) 


where Dyii{z, Mj) represents the probability that the parton labelled I fragments to a photon 
with a momentum fraction z of its own momentum at scale Mj (note that Dy/i{z, Mj) is 
the usual unpolarized fragmentation function, since the hnal state is not polarized). This 
is the non-perturbative fragmentation function which must be obtained from experiment at 
some scale and evolved to Mj using the usual evolution equations. 

In NLO, 0{aa‘^), there are virtual corrections to the LO non-fragmentation processes of 
Eq. (H). as well as the further three-body processes: 


g +q^ g + q + 1 (4.5a) 

g + g ^ q + q + i (4.5b) 

q + q ^ g + g + 1 (4.5c) 

q + g^g + g + 7 (4.5d) 

q +q^q + q + 'y (4.5e) 

q+q^q' + q' + 'J (4.5f) 

q + q' ^ q + q' + (4.5g) 
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In principle the the fragmentation processes of Eq. ( |4.2|) should now be calculated to 
0{al) and convoluted with the NLO photon fragmentation functions whose leading be¬ 
haviour is 0{a/as), but they have not yet been calculated for the polarized case, hence, 
in both the polarized and unpolarized cases, we include the leading order contributions to 
these processes only. 

The direct contribution to the inclusive cross section is given by 




Jvw 1 — u Jvw/v w 




X 


1 dAa. 


0^7 


dv 


5(1 — w) + 


as{iA 
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-AKij^^{s, u, w, Z, M^, Mj) 


(4.6) 


where AKij^y(s,v,w, Z, M^, Mj) represents the higher corrections to the hard subprocess 
cross sections calculated in [O and p is the renormalization scale. 

In this paper we present results for the inclusive prompt photon cross section without 
taking any possible isolation cuts into account. As shown in ref. |^, isolation cuts do not 


have a signihcant effect on the asymmetries, which are the quantities in which we are mainly 
interested here. It is also unlikely that isolation will be necessary at HERA cms energies. 


2. Single Jet Produetion 


We calculate the single jet cross section at LO only. In this case the approximation 
‘parton = jet’ is used and a jet dehnition is not needed. The eight subprocesses listed in Eq. 


( [4.2|) are the same ones contributing to jet production in LO. The differential cross section 
is given by; 


dAa-^ 
' d^pj 




f'" Afl(x„ M^)An(x 2 , 

'vw 1 — U Jvwlv V dv 


-5(1 


w 


( 4 . 7 ) 


where the variables are dehned as before. Recently the polarized jet cross section has 


been calculated in NLO, [Q but for the purposes of this paper we use the leading order 
calculation. 


C. Significance of Results 

We now compare our three sets of polarized distributions, (GGRA, GGRB and GGRG) 
0 generated from the three gluon models discussed in section II, with those of Gehrmann 
and Stirling, [|^] set GSA. Our polarized gluon distributions are smaller than those of 
GSA, which affects the ratio of quark to gluon contributions at small x. This in turn will 
modify the relative contributions to direct-y production and Jet production. (See £g. la). 
The polarized sea peaks at higher x and goes negative at much larger x than the GSA 
distributions, which will further affect the relative contributions at small x. (See hgs. lb 
and Ic). Our valence distributions are also larger, which makes the qq contributions to 
both processes larger. (See fig. Id). For GGRA, the LO and NLO distributions are similar. 
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implying perturbative stability. The result is similar for GGRB and GGRG (figs. Ic and 
Id). 

The comparative results for direct -7 production are shown in hgures 2 through 4. At 
'/S = 200 GeV, the cross section is measurable out to p-p ~ 40 — 50 GeV, given the planned 
luminosities of RHIG (£g. 2a). The asymmetry in direct -7 production differs signih- 
cantly from GSA only at the largest pt values (> 40 GeV), where the PHENIX detector at 
RHIG is at an upper pt limit. However, expected errors in the measurement of the asym¬ 
metry, at PHENIX are much smaller at the lower pt values. Up to pt of about 
15 GeV, the expected error is about ±0.006 and grows to about ±0.080 at 30 GeV. Thus, 
measurement of this asymmetry should be done at all possible pt values, with those at lower 
Pt possibly able to distinguish between the small and large polarized gluon distributions, 
despite their close proximity. (See £g. 2b). At \fS = 200 GeV, the qg processes in GGRA 
dominate (positive AG) much more than the GSA and our other gluon models, which give 
the larger values seen at larger pt- The qq contributions are relatively small at all px 
for all models. See £g. 2c for details. Thus, the qg processes dominate, as expected. In 
£g. 3a, we plot the direct -7 cross section as a function of rapidity, at Px = 15 GeV. The 
cross section remains of reasonable size out to rapidities past | y'^ |~ 1.3. The corresponding 
asymmetries (fig. 3b) are relatively flat out to | |~ 1, after which they exhibit a slight 

rise. They become quite large at rapidities which are out of the range for PHENIX, but 
possibly accessible by the endcap of the STAR detector. It may be advantageous to take 
measurements at | |> 1 at the appropriate ranges of rapidity with both detectors for 

better sensitivity to the asymmetries. 

At HERA energies, a/s ~ 40 GeV, the cross section is reasonable for ranges of px up to 
about 10 GeV (£g. 4a). There are differences in the asymmetry A^^,, which depend upon the 
models of the polarized distributions. However, these differences manifest themselves more 
greatly at px ~ 10 GeV and above. The sizes of A^^. still large enough at the planned px 
values at HERA (about 2 —>■ 8 GeV) to be measured. Since the asymmetries diverge quicker 
at the HERA center-of-mass energy than at higher a/s, there may be a good chance that 
a more precise determination of the polarized gluon can be made. A lot depends upon the 
relative errors expected from such measurements at each p/. (See £g. 4b). Thus, even the 
moderately sized positive AG gives a relatively large asymmetry, Small asymmetries 

would indicate a very small, if not negative AG in this kinematic region. A possible problem 
exists, since the expected experimental errors tend to increase substantially with increasing 
Px- Thus, these two cases may not be distinguishable. It is recommended that be 
measured over a wide range of pr, up to the largest values possible, to test the consistency 
of these models. 

Although the direct -7 cross sections are signihcantly larger at the RHIG energy of a/ 5 = 
500 GeV, the asymmetries are correspondingly smaller and it would be much more difhcult 
to determine the polarized gluon distribution from direct -7 production. Thus, we highly 
recommend that these measurements be done at HERA-iV and at the lower a/ 5 values at 
RHIG. 

Figures 5 through 6 show our predictions for jet production (pp —>• Jet±X) using these 
distributions. The cross sections are measurable for all the planned values of p/. These are 
shown for zero rapidity at y/S = 200 GeV in £g. 5a. The qq and qg processes both dominate 
in our model for A'h^ at large p/. The qg is more dominant at smaller p/, with comparable 




contributions at the limit of the RHIC STAR detector, {pt ~ 25 GeV at y/S = 200 GeV). 
This is consistent with the discussion above. The GSA model has the qg process totally 
dominant with qq small (see fig. 5b). However, their asymmetries are quite a bit smaller than 
ours for > 20 GeV and ours rise with p^, while theirs remain constant. A measurement of 
for different at the smaller yfS region of RHIG may be able to distinguish among these 
models, and give an indication as to which process does indeed dominate in this kinematical 
region. The GS and GGR models are both around the 4 ^ 6 % level at p^ ~ 20 GeV while 
the GGR asymmetries grow with up to the maximum values at RHIG (see fig. 5c). Our 
asymmetries for the other gluon models exhibit the same p^ behavior, except that the values 
are smaller (up to about the 10% level at the limit). Again, it is important that the errors in 
jet measurement at STAR be minimized, in order to distinguish among the models at these 
levels. This can be, at least in part, accomplished with the endcap installed, which will 
provide enhanced jet coverage in this kinematic region. Distinguishing among the various 
models can yield vital information about both the polarized gluons and quark distributions 
in Jet production. (See hgs. 5d and 5e). 

At large \/s ~ 500 GeV, the cross sections are larger than at 200 GeV (see fig. 6 a), but 
the corresponding asymmetries in Jet production are much smaller and cannot distinguish 
between the various models, (fig. 6 b) even at the largest values available. Thus, it is 
recommended that these asymmetries be measured at the lower center-of-mass energies. 


V. DISCUSSION 

We have calculated the cross sections and asymmetries for polarized direct-q production 
and Jet production at RHIG and HERA (for 7 production) cms energies, using the polarized 
distributions generated in reference [Q. Gomparisons of our predictions with those of other 
distributions were discussed in light of the physics that can be extracted from performing 
these experiments. The entire evolution for these predictions was done for LO and NLO 
in x-space using a technique explained in the text and in the Appendix. This gives highly 
reliable results without the use of moments. Both the LO and NLO evolution programs are 
written in FORTRAN and are available from either author via electronic mail. 

We have shown that if the polarized gluon distribution is even moderately positive at 

= 1 GeV, the asymmetry for direct photon production is among the best candidates for 
determining the size of AG. The size of the predicted asymmetries may, however, be less 
relevant if the experimental errors are too large. Jet production is a close contender for 
determining which models of AG give the best agreement with data. Much will depend 
upon the relative errors in the applicable kinematic regions for PHENIX (photon detection) 
and STAR (detection of photons/jets). For details, see the articles by Y. Goto, S. Vigdor, S. 
Heppelman and D. Underwood in ref. At this point, it appears that the projected errors 
for direct -7 production at yfS = 200 GeV may be able to distinguish between predictions 
based upon large and small values of AG. Any sensitivity of polarized direct -7 production 
to kf smearing effects should cancel out when the asymetries are calculated, which should 
eliminate some of the potential theoretical uncertainties. Furthermore, using both PHENIX 
and STAR to measure direct photon production will enlarge the rapidity range, and thus 
increase the sensitivity to extract useful physics information from this process. 
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We have also shown that Jet production provides an excellent mechanism to distin¬ 
guish among polarized gluon models, as well as polarized quark distribution models. The 
corresponding predicted asymmetries are large enough to measure and distinct enough to 
distinguish among many models, including those mentioned in the last section and other 
extreme models, such as a saturated AG model. |]^ It is important to stress that it is nec¬ 


essary to include both triggering and endcap detection in STAR for the fullest coverage of 
X. This will help not only determine the size, but also the shape of AG. The jet production 
rates appear to be sufficient in the appropriate kinematic regions, but detection must be 
complete to extract the maximum physics information. The two processes discussed above 
seem to have greater size and more signihcant differences in the various predicted asym¬ 
metries than production for the same kinematic regions. [jT^ It is thus important that 
experimental errors be minimized for both PHENIX and STAR so that these measurements 
can be made. 

Since HERA-A^ coverage for direct -7 production is at lower cms energy and spans a 
lower pt range, this experiment would be complementary with that of RHIC in the gluon 
kinematical region covered. The sensitivity of HERA-A^ to AG for direct -7 production is 
sufficient to distinguish between the various gluon models. 0 Coupling this with the other 
planned experiments at HERA-A^ would span a sufficient kinematic range to give a very 
good indication as to both the size and shape of AG. ||^ 

There are significant enough differences in the models for the gluon distributions 
to warrant careful study of these processes. Measurement of asymmetries for 


direct photon production and jet production at RHIC and HERA-X are crucial in deter¬ 
mining the size and shape of the polarized gluon distribution. 
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APPENDIX: X-SPACE EVOLUTION 


In this appendix we give details of our x-space evolution of the polarized parton densities. 
We use the dehnitions of ref. for the non-singlet distributions and splitting functions. 
Starting at Qq with Nf = 3 flavours we define non-singlet distributions 

Aqf^ = Aqi ± Ag* (Al) 


and evolve the distributions Aq^ and Aqp' — Aq^ according to the equations 


d{Aqt - Ag+)(x,Q2) 
dings 


^^AP+(x) * (Ag+ - Ag+)(x, Q^) 


(A 2 ) 
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and 


where 


dlng 2 


asm 

2ti 




(A3) 


AF^ = F^ ± AFh. 

qq qq qq 


(A4) 


The singlet distribntion is dehned in Eq. (|3.3|) and the evolntion eqnation is given by 
Eq. (|3.4|). In ref. 0 the dehnition 


AF = AF+ 

^J^qq ‘-^^qq 


AF® 

qq 


(AS) 


is nsed. All the splitting fnnctions are given in that reference. 

For Nf = 3 flavours we have four non-singlet distributions to evolve 

Ts = Am’*' — Ad’*' = Au + Au — Ad — Ad 

Tg = Au'^ + Ad'^ — 2As^ = Au + Au + Ad + Ad — 4As 

Vi = Au — Au 

V2 = Ad- Ad. (A6) 


plus the singlet quark 

S = Am -|- Am -|- Ad -|- Ad -|- 2As 


(A7) 


and gluon, AG, distributions. We have assumed that As = As. The functional forms of 
the input distributions are given in reference [Q. 

We label the coefficient functions for the different sectors by and A^ 

and correspondingly for the B's. In LO the coefficients are given by the relations; 


A(r®i^(a;) = Amix) 


NS,±i 


3C p 4C T. 


Po Po 


ln(l — x) 


+ ^ / ^(1+.-Mrys) - ^ 

Po Jx y Po 


dy yA^^’^iy) - xA^^'^{x) 


y-x 


(AS) 


A^+i(x) = A^{x) 


3Cf 

Po 


4Cf 

Po 


ln(l — x) 


+ 2 /■^{C,.(l + .-)A^(9)-iV^(2z 

Po -Jx y '' 

_ 4C^ F dy yA^{y) - xA^{x) 

Po Jx y y-x 



Al, 


n+l 


4Nc 

Po 


ln(l — x) 


{x) = -A^{x) 1 + 

-1 lAI( 2 -z)A^(y)- — 

Po Jx y [ Po 


ANr 


(2z - l)A°{y) 


iyyAniy) -z:A°{x) 


lx y 


y-x 


(A9) 


(AlO) 
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where z = x/y. 

The NLO coefficients are given by; 


B, 


NS,±, 

n+1 


(X = 


A 


■A 


NS,±, 


47r/3o' 
10 


(X 


B, 


NS,±, 


9 


ln(l — x) — C 


X] 11 + ^(3 + 41n(l - x) 
Po 

TT^ \ /17 

- — + 6C - iVj — + 


+ 


C, 


TT/^O 

-3C + 


TT 


2Cf 


Jx y 


J “ V24 18 

4CFyB:^^’^y)-xB^^’^x) 


Nf\ — + — 
^\12 9 

67 

Y 


Po 


(1 + ^ 


Po 
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Figure Captions 


[1] (a) Polarized gluon distribution xAG{x, Q"^) at = 100 GeV^ showing the 
GSA, GGRA, GGRB and GGRG parametrizations. (b) same as (a) for the Au- 
quark distribution, (c) Same as (a) and (b) for the strange quark distribution. 
The dotted curve near the solid (GGRA) is for the LO version of GGRA. Same 
as (c) for the valence u-quark distribution. Again the dotted curve near the solid 
curve is the LO version of GGRA. The valence distributions for GGRA, GGRB 
and GGRG cannot be distinguished. 

[2] (a) pT-distribution for unpolarized prompt photon production ed. y'^ = 0 for 
y/S = 200 GeV. (b) Asymmetry Alj^ for the cross section in (a) as predicted 
by the GSA, GGRA, GGRB and GGRG polarized parton distributions, (c) 
Asymmetries from (b) for GSA, GGRA amd GGRG showing the contributions 
from the qg and qq initial state contributions. 

[3] (a) Rapidity distribution for unpolarized prompt photon production at = 15 
GeV and y/S = 200 GeV. (b) Asymmetries for the rapidity distribution as 
predicted by various polarized parton parametrizations. 

[4] (a) pr-distribution for unpolarized prompt photon production at = 0 for 

y/S = 40 GeV. (b) Asymmetry ^he cross section in (a) as predicted by 

the GSA, GGRA, GGRB and GGRG polarized parton distributions. 

[5] (a) pT-distribution for the unpolarized jet cross section at = 0 for y/S = 200 
GeV. (b) Asymmetry A'h^ for the total and for various initial state contributions 
as predicted by the GSA distributions, (c), (d) and (e): Same as (b) for the 
GGRA, GGRB and GGRG parametrizations respectively. 

[6] (a) pT-distribution for the unpolarized jet cross section at p*^ = 0 for y/S = 500 
GeV. (b) Asymmetry A'lj^ for the cross section as predicted by the GSA, GGRA, 
GGRB and GGRG polarized distributions. 
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